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Abstract. We prove Strichartz estimates for the absolutely continuous 
evolution of a Schrodinger operator H — {iV + ^4)^ + V in R", n > 
3. Both the magnetic and electric potentials are time-independent and 
satisfy pointwise polynomial decay bounds. The vector potential A{x) is 
assumed to be continuous but need not possess any Sobolev regularity. 
This work is a refinement of previous methods, which required extra 
conditions on divy4 or |V| 2 A in order to place the first order part of the 
perturbation within a suitable class of pseudo-differential operators. 



1. Background 

This paper continues an investigation into the dispersive properties of 
Schrodinger operators taking the form 

H = {i^7 + A{x)f + V{x) = -A + i{A-V + V ■ A) + \A\^ + V 

=A + L{x) 

where A{x) is a vector field in M", n > 3, and V{x) is a scalar function. The 
associated evolution equation 



(2) 



iut{t, x) = Hu{t, x) 
u{0, x) = uq{x) 



models the motion of a single charged particle within an ambient electro- 
magnetic field, with V and A serving respectively as the electrostatic and 
magnetic potentials. Such operators also appear routinely when a nonlinear 
Schrodinger equation is linearized around a nonzero static solution. 

The starting point for estimates in the free case {L = 0) is an explicit 
formula for the propagator e**^ based on Fourier inversion. 

n(t,x) = (47rit)""/^ / e'^'^^ uo{x) dx 
It is immediately clear that u = e**^iio satisfies a family of dispersive bounds 

(3) Mt, ■)\\q<\tr"^'^'~'^ho\\g', 2<(?<00, 
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the p = 2 case following from Plancherel's identity and the p = oo case from 
evaluation of the convolution integral above. Each of the dispersive bounds 
with 2 < q < can be incorporated into a TT* argument to prove one of 
the Strichartz inequalities, 

(4) ||e'*''^^ollL?L| < holb, | + f = t, 2<p,g<cx). 

This line of argument requires nothing more than the Hardy-Littlewood- 
Sobolev inequality (except for the p = 2 endpoint case which is more deli- 
cate |14]). and the resulting estimate have played a central and the resulting 
estimates have played a central role in the well-posedness theory of nonlinear 
Schrodinger equations from their inception. It would therefore be desirable 
to prove an analogous statement for solutions of ([2]) with more general H. 
Unfortunately, it may not suffice to replace the Laplacian in ([2]) with another 
Schrodinger operator H because of the possible existence of point spectrum. 
Each eigenvalue of H gives rise to a solution of the form n(t, x) = e~^'^^'ip{x), 
negating the bound 

for every p < oo. 

For short range self-adjoint perturbations of the Laplacian, the spectrum 
of H should remain absolutely continuous along the positive real axis, with 
discrete negative eigenvalues possibly accumulating at zero. No eigenvalues 
are embedded within the continuous spectrum flS]. In many cases it is 
possible to prove that the linear Schrodinger evolution decomposes into a 
discrete sum of bound states, plus a radiation term that enjoys the same 
dispersive properties as a free wave. Our goal is to expand the class of 
potentials for which this behavior is known to occur. 

Following the notation in Chapter XIV of [lOj . define the function space 

oo 

(5) \\f\\B:=J2'^'\\f\\L^iD,), ll/b* :=sup2-i||/||i2(^^.), 

where Dj = {x £ : \x\ ~ 2-'} is a decomposition of into dyadic 
shells for j > 1 and Dq = {\x\ £ M" : |x| < 1}. These share numerous 
characteristics with weighted L^, and we designate Bs to be the associated 
homogeneous Sobolev spaces with norm 

(6) ll/b. :=ll|Vr/lb, ll/lk:=ll|VlVb- 

One other Banach space related to B will also play an important role during 
the discussion. Let 

°° 1/2 

(7) lblU:=(E2'^'Nli-(D.)) 

j=0 

Observe that pointwise multiplication by a function in X maps B* to the 
weighted space {x)~^ L'^(W^), and also maps {x)!? to B. 
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Theorem 1. Let H = (iV + A)"^ + V he a magnetic Schrddinger operator 
on M", n > 3, whose scalar and magnetic potentials are bounded functions 
satisfy the conditions 

(CI) {xfV e L'^iM''), and limsup |x|2|y(x)| = 

(C2) A is continuous, and \\A\\x < oo 

Let Pac{H) represent the orthogonal projection onto the absolutely continu- 
ous spectrum of H . If zero is not an eigenvalue or a resonance of H , then 
the Strichartz inequalities 

(8) \\e-^'^PaciH)uoh^^Ll < hoh, f + f = t 

are valid for each exponent p > 2, as well as the Kato smoothing bound 

(9) Ux)\V\^l\-''" PUH)uo\y^L, < Ib'llill^^olb. 

It is not possible to follow precisely in the footsteps of the free case. 
In particular, no replacement for ([3]) is known in the case A ^ 0, and to 
achieve the full range of exponents likely requires additional regularity as- 
sumptions [9], [3]. There are two alternative approaches to this problem, 
both of which attempt to confine perturbative arguments to the more for- 
giving L^ setting. 

The first method is rooted in pseudo-differential calculus, decomposing 
the solution into wave packets and finding suitable parametrices to de- 
scribe their respective trajectories. This techniques were applied to mag- 
netic Schrodinger operators as early as [17]. More recent advances have 
highlighted the flexibility to work on manifolds other than R" (as in. [2]) 
and with time- varying coefficients [l6j . 

We instead adopt the framework outlined in [18], where Strichartz esti- 
mates follow directly from a pair of Kato smoothing bounds, one each for 
the free and perturbed propagators. The starting point is to express the 
perturbed evolution according to Duhamel's formula, with H = —A + L. 

(10) e-'^^uo = e'^^uo - i f e'^^''^^ Le~''" uq ds 

Jo 

The free evolution term is controlled by ([1]). To prove a Strichartz estimate 
for the integral term, it typically suffices to find a factorization 

J 

i=i 

for which Zj is a smooth perturbation relative to H on its absolutely con- 
tinuous subspace and Yj is smooth relative to the Laplacian. Following 
Kato's pL3j theory of smoothing, this reduces to proving uniform estimates 
for the resolvents of both H and —A. 

It is well known that Yj must not have order greater than ^, and one 
expects the same to be true of Zj. The natural factorization of L = ^ • V 
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would therefore distribute half of a derivative to each of Yj and Zj. Some 
regularity of the coefficients A(x) is needed in order to commute half a 
derivative across them. 

To handle less smooth magnetic potentials, we take advantage of the par- 
abolic nature of the Schrodinger equation to replace powers of the gradient 
with derivatives in the time direction. Since A{x) is time-independent, it 
commutes with operators such as {dt) without apparent difficulty. There is 
a price to be paid later for this convenience: some of the precise structure of 
Duhamel's formula (i.e. < s < t) is lost in the process. Another apparent 
sacrifice is that the polynomial decay condition ()C2p falls one-half a power 
short of the natural scaling law A{x) ~ 

In the next section we choose a factorization for the first-order perturba- 
tion found in ([T]) and prove preliminary smoothing bounds for each term. 
The task is not complicated, as many of the underlying resolvent estimates 
appear in the literature in precisely the form needed. Section [3] contains the 
direct proof of Theorem [TJ The key ingredient is a set of modified smoothing 
estimates whose form is dictated by the domain of integration in (jlOp . The 
technical issues raised by our use of fractional time-derivatives are addressed 
here. 



2. Resolvent and Smoothing Estimates 

The magnetic and scalar potentials collectively perturb the Laplacian with 
a differential operator of the form 

L = i{A -V + V ■ A) + \A\'^ + V 

where our factorization of choice is to take 

Yi = {x)-\dt)-*, Y2 = i{x){dtr\A-V, Y^ = \V +\A\^\^'\ 
Zi = Y2, Z2 = Fi, Zs = sign(y + \A\^)Y^. 

Operators Yi and Y2 produce vector-valued functions, and could be further 
split (if necessary) into their three coordinate directions. The first set of 
smoothing estimates for Yj relative to the Laplacian are now easy to obtain. 

Lemma 2. Suppose V{x) and A{x) are hounded and satisfy (jCip . ()C2p . 
Each of the operators Yj, j = 1,2,3, satisfies the bound 

(12) r,e**^/llL2(MxR") < II/II2 
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for all functions f € L^(M"). To prevent ambiguity, the first two of these 
mapping estimates are interpreted as follows. 

(13) / {x)-'\\e^'''f{-,x)\\\ dx<\\fg 



with the intermediate function e**^/ taking values at all times t G M, both 
positive and negative. 

Proof. Let A represent the Fourier variable dual to t. After setting up the 
TT* operator and taking the Fourier transform in T, the smoothing bounds 
for Yi, Y2, I3 are essentially equivalent to these properties of the free resol- 
vent. 

(i?o+(A) - i?o (A)) {xr'fW, < (A)-^ II/II2 
(15) 1 

\\{x)gV{R+iX) - R^iX))V{x)gf\\^< {\)^\\gfM\2 

uniformly over A > with any exponent a > 1. We have adopted the 
shorthand notation for resolvents 

Rt{X)=lim{H-{X±ie))-' 

(16) 

R^iX)=]im{-A-{X±ie)) 

The first inequality can be found in [1], and a particularly sharp version 
appears in [T9]. The second is a standard result in [10], together with the 
observation that Ai?^(A) = -/ - Ai?^(A). □ 

In order to show that each Zj is a smooth perturbation relative to (the 
absolutely continuous part of) H, it suffices to verify that the resolvent 
estimates in (fT5|) continue to hold when -R^(A) is replaced by R^{X). It is 
convenient to break the problem into three separate regimes according to 
whether A<Cl,A~l,orA3>l. The latter two cases have been considered 
at length elsewhere, so we are content to collect these results into a single 
statement. 

Lemma 3. Let A be a vector field and V a scalar function satisfying con- 
ditions ()Cip - (|C2p . Given any number Aq > 0, there exists a constant 
C{L,Xo) < 00 such that 

^^^^ \\{xy'{Rt{X) - Rj^{X)){x)-'f\\^<CiL,Xo){X)--2\\fh 

\\{x)gV{Rl{X)-Rl{X))V{x)gf\\^<C{L,Xo){X)'2\\g\\j, 



uniformly over all A > Aq- 
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Proof. According to Proposition 4.3 of [7], there exists a number Ai(L) < oo 
such that the perturbed resolvents R^{\) satisfy the bounds 

\\Ri{X)\\B^B' <CX-^ 
||i?^(A)|b_,^B. <CA5 

for ah A > Ai. A uniform bound over the interval A € [Ao,Ai] is proved 
as part of Theorem 1.3 in [11]. It is assumed there that no eigenvalues 
are embedded in the positive half-line, a condition which was subsequently 
shown in [15j to hold for the entire class of potentials under consideration. 

□ 

Remark 4. The mapping bound from B to B* is considerably stronger than 
what is needed to prove Lemma El Alternatively, by projecting onto the 
range A € [Ao,oo) rather than the entire continuous spectrum, one can 
obtain smoothing estimates while assuming less decay of the potentials. It 
should suffice to let A, V belong to the space 

oo 

X := e L°° : lbllL-(D,) < oo}- 
j=0 

The remaining short interval A G [0, Ai] is handled by a compactness 
argument. By applying resolvent identities, the difference i?^(A) — RJ^{X) 
can be expressed as 

Rt{X)-Rl{X) = {I + R^{X)L)-'{R+{X)-R^{X)){I + LR^{X))-' 

It is already established (see (llSp ) that the difference of free resolvents maps 
the space Y := {x)-^L^{W) + equipped with the norm 

(18) ||/||y = inf{||(x)/i||2 + ||/2b_, : /i + /2 = /} 

to its dual Y* = {x)L'^ fl B^ with uniformly bounded operator norm. We 
therefore need only to show that (/ + Rq {X)L)~^ exists as a uniformly 
bounded family of operators on Y* over this range of A. The same will be 
true of the other inverse by duality. 

Lemma 5. Let V and A satisfy conditions ()Cip - ()C2p and suppose that the 
Schrodinger operator H does not have an eigenvalue or resonance at zero. 
Then 

sup II (/ + RQ{X)Ly'^\\y,^y, < OO. 

AG[0,Ao] 

Proof. Any potential V € satisfying (IClh or A £ X can be approx- 

imated in the appropriate norm by a function with compact support. As a 
consequence, L acts as a compact operator taking Y* to Y. Meanwhile, the 
resolvent Rq^{X) maps Y back to Y* . The end result is that / -|- Rq{X)L is 
always a compact perturbation of the identity. 
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The Fredholm Alternative asserts that (/ + Rq{X)L) must possess a 
bounded inverse unless there there is a nontrivial null-space consisting of 
functions / G Y* that solve the equation / = — i?(J"(A)L/. Any such / 
would also be a distributional solution of the eigenfunction equation 

-Af + Lf = Xf. 

The combined efforts of [15] and [11] rule out their existence for each A > 0. 
To avoid them when A = we must make include an a priori assumption 
that zero is not an eigenvalue or resonance of H. The need to exclude 
resonances stems from the fact that Y* includes functions that do not decay 
rapidly enough to belong to L^{W^). 

Remark 6. In [12] the presence of an eigenvalue at A = is shown to nullify 
some Strichartz estimates even after projecting away from the associated 
eigenvector with Pac{H). 

Now that (/ + Rq {X)L)~^ exists at each A G [0, Aq], the uniform bound 
on their norms follows from continuity with respect to A. More precisely, 
the resolvents Rq{\) enjoy the weak-* continuity property described below. 
Given two functions f,g€Y and A E [0, oo), 

hm{R+{X)f,g) = {R+{A)f,g). 

The main idea here is again that / and g can be approximated by compactly 
supported functions, and that the integration kernel of Rq{A) — Rq{X) van- 
ishes pointwise on bounded sets. 

Now suppose there were a sequence of values A^ for which the norm of 
(/ -|- Rq {Xn)L)~^ grew without bound. This would imply the existence of 
unit functions /„ € Y* for which ||/n + Rq {Xn)fn\\Y* — > 0. By passing to 
a subsequence we may assume that A„ ^ A and that /„ converges in the 
weak-* topology to a function /oo € Y*. 

Recall that L is a compact operator, which means that L/„ should con- 
verge to Lfao € y in norm. By the weak-*continuity of resolvents, it follows 
that Rq {A)Lfoo is the weak-* limit of Rq {Xn)Lfn- This sequence is suffi- 
ciently close to — /„ that — i?^(A)L/oo must be a weak-* limit of /„ This 
makes /oo a solution to the eigenfunction equation, so /oo = 0. 

The contradiction arises because now we have |[L/„||y — > 0, which would 
cause ll/n + Ro{Xn)fn\\Y* 1; violating the part of the construction where 
||/„ + i?+(A„)L/„|| ^0. □ 

These resolvent estimates can now be combined to prove the second Kato 
smoothing estimate for our factorization of L. 

Lemma 7. Suppose V{x) and A{x) are bounded and satisfy (jCip . ()C2p . 

Each of the operators Zj, j = 1,2,3, satisfies the bound 

(19) ||Z,e-**^P,e(^)/||L2(KxRn) < ll/lb 
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for all functions f € L^(M"). To prevent ambiguity, the first two of these 
mapping estimates are interpreted as follows. 

£(x)2||^-V(e-*^P,e(i/)/)(-,x)||^_i(^^dx<||A||2,||/||2 

with the intermediate function e^^^^ PaciH)f taking values at all times t £ 
M, both positive and negative. 

The Kato smoothing conclusion ([9|) in Theorem [1] is a direct interpolation 
of the two inequalities above. 

3. Proof of Theorem [T] 

The intended purpose of Lemmas [2] and [7] was to control the integral 
term in Duhamel's formula (llOp . For each term in the factorization of L, 
one would like to apply (119p to obtain a function in L^(M x M"), followed 
by the dual form of (|12p to return it to L'^(M"), with the free Strichartz 
estimate (j4]) then giving a map into L^L%. 

The domain of integration in (llOp is clearly stated as < s < t. Unfortu- 
nately, the interpretation of fractional derivatives and time-propagation in 
both (|12p and (|19p is just as clearly two-sided in s. The desired Strichartz 
estimate dH) therefore follows from two modified smoothing bounds which 
include a sharp cutoff to preserve the triangular region of integration. 

Lemma 8. Suppose V{x) and A{x) are bounded and satisfy (IClj) . (IC2p . 
Each of the operators Zj, j = 1,2,3, satisfies the bound 

(20) ||Z,(l{,>o}e-^)P,e(i/)/||^,^, < II/II2 

for all functions f G L^(M"). 

Lemma 9. Under the same hypotheses, each of the operators Yj, j = 1,2,3, 
satisfies a bound 



(21) / (l{i>,}e^(*-^)^)i-*5(s,x)ds 

Jm. 



for each Strichartz pair | + ^ = fi P>2, and every 17 G L^(M x M"). 

Proof of Lemma O For Z3 this is a trivial consequence of (jl9p . For both Zi 
and Z2 it is also an immediate consequence, as pointwise multiplication by 
l{s>o} is a bounded operator on every Sobolev space i/", |a| < i. □ 

Proof of Lemma\^ For the pointwise multiplication operator 13, the argu- 
ment in [1^ requires no modification. The dual form of (jl2p combines 
with (jl]) to yield a bound 



/ e'^'-''>''Y^g{s,x)ds 
Jm. 
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for all Strichartz pairs {p,q) including the endpoint (2,^^). The Christ- 
Kiselev lemma [5] (in particular the version appearing in [20]) confirms that 
the same bound is true if the integral is only taken over the triangular region 
t > s, provided there is a strict inequality p > 2 in the exponents. 

When the same logic is applied to Yi , the first integral estimate takes the 
form 



{{d,)le^('''^^){x)-'g{s,x)ds 
Directly applying the Christ-Kiselev lemma leads to the further bound 

The statement of (f2T|) for Yi instead concerns the operator inequahty 



[ {ds)Hl{t>s}e'^'-'^^){x)-'gis,x)ds 

JR 



SO we are left to bound the difference between the two. This consists of two 
terms, 



IS 



- l{i>,}(a,)i (l{i<,}e^(*-^)^)) {x)-^g{s, x) ds. 

It suffices to control the first term, as the second will behave identically by 
symmetry. After the initial multiplication by {x)~^, the rest of the operator 
is a convolution in M x M" with a kernel Ki[t,x) satisfying the bounds 



(22) 



\Ki{t,x)\ 



-{t\x\ 



g-t |a;|-"+2(2;)- 



if i < 1, 
if t > 1. 



A direct computation then shows that Ki belongs comfortably to both 
L^Lx^'^ and L1l\, thus convolution against Ki maps L2(M X W-) (even 
without the extra weight (x)~^) to any of the mixed-norm Strichartz spaces. 

The analysis for Y2 is quite similar. Once again there is a bulk term that 
can be controlled by combining the dual form of (jl2p with ^ and the Christ- 
Kiselev lemma. There is a nonzero remainder generated by commuting 
multiplication by the cutoff lt>s with fractional integration in s. In this 
case its precise form is 



-l|t>.}(a,)-^(l{i<,}Ve^(*-^)^ 



A{x){x)g{s, x) ds. 
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As before, multiplication by A{x) ^ 1 is composed with a convolution oper- 
ator whose kernel K2{t,x) satisfies the bounds 



(23) 



\K2{t,x)\ 



< 



{t\x\ 



{x)- 



if i < 1, 
if i > 1. 



The kernel K2 belongs to both L^L^. and L|L"/^'°°, which is still sufficient 
to generate a bounded map from L^. to each of the Strichartz spaces. 

The only task remaining is to verify the asserted size estimates (|22p 
and ()23p . Each one involves a fractional derivative or integral, so they are 
best handled via the Fourier transform. Recall that the convolution kernel 
for the Schrodinger propagator can be expressed as 



L{t>0}« 



itA 



(x) 



I 

Jm. 



-itX 



(A){x) dX 



therefore the kernel Ki (t, x) is defined to be 



(24) 



Ki(t,x) = l|i<o} [ {X)h-''^R-iX)ix)dX. 

JR 



For any t < 0, the free resolvent R (A) has an analytic continuation in A to 
the lower halfplane. Moreover, its kernel has the asymptotic bound 



1 7i-3+2|q 



\D''R-{X){x)\ < |2;r"+2-l°lel^'l^"^(^'^')(|x|Ai) 

for derivatives of order |a| = 0, 1. 

The smooth power function (A) 4 is holomorphic in the halfplane with a 
segment removed extending from —i downward along the imaginary axis. 
This allows the contour of integration to be moved from the real axis to 
the two sides of this imaginary segment. To be precise, the original Fourier 
transform existed only in the distributional sense, so we must first mollify 
the behavior at infinity with a function such as ^^^^ and take limits as 
k 00. 

On this new path, change variables so that A = —ifi^. The result is that 

/oo 
nil - ^x^)^e~*^''R-{-i^?)(x) dn 

/oo _ 
^2e-*^'e-'^N'^(|x|/i)Vd/i 

/oo ^ 
nil - n'^)-se-*''^VR-{-in^){x) dfi 

POD 

-"+1 ' /i3e-*^'e-^l^-|'^(|x|//)'^(i/i 
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where c = ■\/l/2 is a fixed constant. One possible bound comes from domi- 
nating the factor e~'^l^l'^(|a;|//)^ by a constant, which leaves 

/oo 
/x5e-*^' dii < min {t'^ , e"*) |xr"+2 

/oo 
/X2 e-^"' dn < min (t" 4 , e"*) |xr"+i 

These appear to be optimal when |xp < t < 1, and are adequate for our 
purposes when < 1 < t. Another alternative is to dominate e"*'' by e~*, 
in which case 

/oo _ 

/oo _ 

These bounds appear to be sharp when t < |a;p < 1, but are also sufficient 
to complete the argument whenever t,\x\ > 1. 

□ 
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